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Abstract 



a^ 

We discuss the recently achieved Bose-Einstein condensation for alkah-metal atoms in magnetic 
pj . traps. The theoretically predicted low-energy collective oscillations of the condensate have been 
^ ! experimentally confirmed by laser imaging techniques. We show by using Poincare sections that 

_j \ at higher energies non-linear effects appear and oscillations become chaotic. 
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cn ' 1 Introduction 

o 

Q^ ■ Finite many-body quantum systems, both fermionic and bosonic, exhibit collective oscillations, 
G^ . which are nowadays experimentally detected with sophisticated devices. In many cases, these 
collective oscillations can become strongly aperiodic and also chaotic [1-3]. 

^: 

Y ■ III this short report, we analyze the collective modes of a trapped Bose condensate of alkali 
O \ atoms. The Bose-Einstein condensation (BEC) is the macroscopic occupation of the ground- 
r-| \ state of the system of bosons. From 1995 we have experimental results interpreted as an 
J-> ■ evidence of BEC in dilute vapors of confined alkali-metal atoms (^^Rb, ^^Na and ^Li) [4-6]. The 
experiments with alkali-metal atoms generally consist of a laser cooling and confinement in an 
^ , external potential (a magnetic or magneto-optical trap) and an evaporative cooling (tempera- 
ture of the order of 100 nK) [4-6] . Nowadays a dozen of experimental groups have achieved the 
BEC by using different geometries of the confining trap and atomic species. 



The dynamics of the Bose condensate can be accurately described by the Gross-Pitaevskii 
equation [7] of mean-field approximation. The theoretically predicted low-energy collective 
oscillations of the condensate have been experimentally confirmed by laser imaging techniques 
[8]. In this short report we show that at higher energies non-linear effects appear and eventually 
the collective oscillations become chaotic. 



2 Mean-field theory of BEC 

The many-body problem of N identical bosonic atoms in an external trapping potential can be 
formulated within the non-relativistic quantum field theory [9-11]. The Lagrangian density of 
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the Schrodinger field is given by 

£ = #.+ (r,«) 



'4 + £^^ - ^w 



V^(r,t) 



--y rfV V^+(r,t)V^+(r',t)\/(|r-r'|)V;(r',t)^(r,t) , (1) 

where U{r) is the external potential and V^(|r — r'|) is the interatomic potential. The quantum 
bosonic field ip{r) satisfies the standard equal-time commutation rules. 



The Lagrangian is invariant under the global f/(l) gauge transformation 

V^(r,t)-.e^>(r,t), 
which implies the conservation of the total number of particles 

N = [ £r V^+(r,t)7/'(r,t) . 



(2) 



(3) 



The Bose-Einstein condensation (BEC) is the macroscopic occupation of the A^-body ground- 
state |0 > of the system. To study BEC a useful mean-field prescription is to separate out the 
condensate contribution to the bosonic field operator in the following way 



^(r,t)=0(r,t) + S(r,t), 



(4) 



where 0(r,t) =< 0|-?/'(r, t)|0 > is the so-called macroscopic wavefunction (or order parameter) 
of the condensate, and S(r,t) is the fiuctuation operator, such that S(r, t)|0 >= 0. Note that 
this prescription breaks the U{1) global gauge symmetry of the system [9-11]. 

Alkali vapors are quite dilute and at zero temperature the atoms are practically all in the 
condensate [4-6]. Thus we can neglect the quantum depletion due to the the operator S and 
the macroscopic wavefunction is normalized to the total number of atoms. Moreover, the 
range of the atom-atom interaction V{r) is believed to be short in comparison to the typical 
length scale of variations of atomic wave functions. The atom-atom interaction is usually 
replaced by an effective zero-range pseudo-potential, V{r) = g5^{r), where g = Airh ag/m 
is the scattering amplitude and a^ is the s-wave scattering length. This scattering length 
is positive (repulsive interaction) for ®^Rb and ^^Na but negative (attractive interaction) for 
^Li. Within these approximations, the Lagrangian density becomes a local function of the 
condensate wavefunction, namely 



£ = 0*(r,t) 



'4t + li^^ - ^(^) 



(r,t)--g\<Pir,t)\' 



(5) 



By imposing the least action principle one obtains the following Euler-Lagrange equation 



. (9 , , , 



^V2 + f/(r)+^|0(r,t)p 
2m 



(r,t) 



(6) 



which is a nonlinear Schrodinger equation and is called time-dependent Gross-Pitaevskii (GP) 
equation [7] . Note that the GP equation is nothing but the mean-field (Hartree) approximation 
of the exact time-dependent Schrodinger equation of the N-body problem, where the totally 
symmetric many-particle wavefunction \1' of the system is decomposed in the following way 

^(ri,r2...,rAr,t) = 0(ri, t)0(r2, t)...0(rAr, t). 



3 Hydro dynamic equations of the Bose condensate 

The complex macroscopic wavefunction 0(r, t) of the condensate can be written in terms of a 
modulus and a phase, as follows 



0(r,t) = v/p(r,t)e^^^^'*^ (7) 

The phase S fixes the velocity field v = {h/m)'VS. The GP equation can hence be rewritten 
in the form of two coupled hydrodynamic equations [12] for the density and the velocity field 

^p + V-(vp)=0 (8) 

^-.^v{L,^,,.—-V'Vr'+—)=0. (9) 

If the repulsive interaction among atoms is strong enough, then the density profiles become 
smooth and one can safely neglect the kinetic pressure term in fi^ in the equation for the 
velocity field, which then takes the form 

m^^^r + v(u + gp+^]=0. (10) 

In the current experiments with alkali metal-atoms, the external trap is well approximated by 
a harmonic potential 

Uir) = -iiu^x' + U2y' + Usz'). (11) 

The ground-state solution (v = 0) of Eq. (10) is given by 

piv)=g''[f^-Uiv)] , (12) 

in the region where /i > U{r), and p = outside. The normalization condition on p(r) provides 
p = {huJh/2)(\bNas/ ahY^^ , where ojh = {ujiuj2UJ^Y/^ and ah = {h/muJhY/'^. 

An analytic class of solutions [13] of the time-dependent hydrodynamic equations (8) and (10), 
which are valid when the condition Nag/ ah >> 1 is satisfied, is found by writing the density 
in the form 

p(r, t) = ao{t) — ai{t)x'^ — a2{t)y^ — a^(t)z^ (13) 

in the region where p(r, t) is positive, and the velocity field as 

v(r,t) = h^[h{t)x^ + h{t)y^ + h{t)z^] . (14) 

The coefficient a^ is fixed by the normalization of the density ao = {lhN/S'Ky'f^{axayazY/^. 
By inserting these expressions in the hydrodynamic equations one finds 6 coupled differential 
equations for the time-dependent parameters ai{t) and hi{t). By introducing new variables g,, 
defined by ai = muf {2gqfqiq2q3)~'^ , the hydrodynamic equations give Oj = (ji/qi and 

qi + u;^q, = i — , (15) 

qiqiq2q3 

with i = 1,2,3. The second and third terms give the effect of the external trap and of the 
interatomic forces, respectively. It is important to observe that, using the new variables g^, 
the equations of motion do not depend on the value of the coupling constant g. In terms of 
qi the mean square radii of the condensate are < xf >= {2fi/muJi)qf and the velocities are 

<vf>=i2p/mu;,)qf [13]. 



4 BEC Collective Modes and Chaos 

The three differential equations (15) are the classical equations of motion of a system with 
coordinates qi and Lagrangian given by 

L = ^(^r'g? + ^2^ql + ^fql) -7^{ql + qi + ql) • (16) 

2 2 gigags 

This Lagrangian describes collective modes of the Bose condensate for Nag/ah » 1 [13]. As 
stressed previously, in such a case the collective dynamics of the condensate does not depend 
on the number of atoms or the scattering length. The minimum of the effective potential is at 
Qi = 1, i = 1, 2, 3. The mass matrix M of the kinetic energy and the Hessian matrix A of the 
potential energy at the equilibrium point are given by 





M = \ LU2 \ and A = 1 3 1 . (17) 



The low-energy collective excitations of the condensate are the small oscillations of variables 
gj's around the equilibrium point. The calculation of the normal mode frequencies fl for the 
motion of the condensate is reduced to the eigenvalue problem A — Q'^M = 0, which gives 

n^ -3 {cuf + tui + cul) n^ + 8 {uf ojI + ojI ojI + ul ojfj n^ - 20 uf cul uj = . (18) 

Note that this formula has been recently obtained by using a variational approach [14] and also 
by studying hydrodynamic density fluctuations of the condensate [15]. For an axially symmetric 
trap, where ui = uj2 = uj±, the previous equation gives 

fii,2= (2 + ^A2±i(l6 + 9A^-16A2)'^') uj^ , n^ = V2 u^ , (19) 

where A = oo^/uo^ is the asymmetry parameter of the trap [12]. Observe that the experimental 
results obtained on sodium vapors at MIT (A = \/S) are in good agreement with the theoretical 
values predicted by (19) [16]. In the case of an isotropic harmonic trap (A = 1) with frequency 
cu, one obtains VLi^2 = a/5c<j, Vl^ = ^/2 to [12]. 

In most experiments the confining trap has axial symmetry [4-6]. Let us analyze this case in 
detail. Because of the axial symmetry, we can impose qi = q2 = q±- Moreover, by using the 
adimensional time r = ui±t, the Hamiltonian of the BEC collective modes becomes 

H = pI + hy, + ql + l-ql + -^ , (20) 

^ ^ q±q3 

where pj_ = dq±/dT and p^ = \~'^dq^/dT are the conjugate momenta. Note that the condition 
gi = g2 restricts collective modes to monopole oscillations, where the third component of the 
angular momentum, that is a good quantum number, is zero [12]. Near the minimum of the 
potential the trajectories in the phase-space are periodic or quasi-periodic. On the contrary, 
far from the minimum, the effect of nonlinearity becomes important. As the KAM theorem 
[17] predicts, parts of phase space become filled with chaotic orbits, while in other parts the 
toroidal surfaces of the integrable system are deformed but not destroyed. 

We use a symplectic Euler method (leap-Frog) to numerically compute the trajectories. The 
time-step is Ar = 10~^ and the energy is conserved to the sixth digit. The conservation 



of energy restricts any trajectory of the four-dimensional phase space to a three-dimensional 
energy shell. At a particular energy, the restriction q_L = I defines a two-dimensional surface 
in the phase space, which is called Poincare section. Each time a particular trajectory passes 
through the surface a point is plotted at the position of intersection (^3,^3). We employ a 
first-order interpolation process to reduce inaccuracies due to the use of a finite step length. 
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Figure 1: Poincare sections with A = \/2. Each panel is at a fixed energy. From left to right 
and from top to bottom: x = 1%! X = 12%, X = 28%, x = 36%. x is the relative increase of 
the energy with respect to the ground-state (minimum of the potential energy). 

In Figure 1 we plot Poincare sections of the system with A = \^. In each panel there is a 
Poincare section with a fixed value of the energy of the system. At each energy value, we have 
chosen different initial conditions [q±{0),q3{0),p±{0),p3{0)] for the dynamics. Actually j9_l(0) has 
been fixed by the conservation of energy. Integration time is 400 in adimensional units, that is 
less than 1 second (the life-time of the condensate is about 10 seconds). Note that the CPU time 
to calculate a Poincare section with a dozen of initial conditions is about 30 seconds. Chaotic 
regions on the Poincare section are characterized by a set of randomly distributed points and 
regular regions by dotted or solid curves. Let x be the relative increase of the energy with 
respect to the ground-state (minimum of the potential energy). For x = 1% and x = 12% 



the trajectories are still all regular but for x = 28% there is chaotic sea. For x = 36% most 
trajectories are chaotic. 

It is important to observe that a strong enhancement of nonlinear effects and eventually chaos 
can be obtained not only by increasing the energy but also by changing the anisotropy A of the 
trap. In fact, as shown in Ref. [13], for special values of A, frequencies of different modes, or of 
their harmonics, can coincide. We plan to investigate in detail the onset of chaos for different 
configurations of the external trap. 



5 Conclusions 



In this short paper we have discussed the mean-field equations which describe the collective 
motion of a trapped weakly-interacting Bose condensate. We have shown by using Poincare 
sections that for large energy values the system becomes chaotic. An important question is 
the following: Can BEC chaotic dynamics be experimentally detected? In our opinion the 
answer is positive. Nowadays non-destructive images of the dynamics of the condensate can 
be obtained. The radius of the condensate as a function of time can be detected and its power 
spectrum analyzed. In fact, one finds the the power spectrum of a chaotic signal is much more 
complex than for a regular one. Typically, one see few peaks in the regular signal and many 
peaks surrounded by a lot of noise in a chaotic one. Finally, we observe that various initial 
conditions for the collective dynamics of the condensate can be obtained by using laser beams 
or by modulating for a short period the magnetic fields which confine the condensate. 
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